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Abstract 

We consider a modification of gravity at large distances in a Brane Universe which was 
discussed recently ]I| - In these models the modification of gravity at large distances 
is ultimately connected to existence of negative tension brane(s) and exponentially small 

tunneling factor. We discuss a general model which interpolates between Bi-gravity H h 

model |jj and GRS model f2j. We also discuss the possible mechanism of stabilization of 
negative tension branes in AdS background. Finally we show that extra degrees of freedom 
of massive gravitons do not lead to disastrous contradiction with General Relativity if the 
stabilization condition / dy V—G^ (T^ — 2T|) = || is implemented. 
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Recently, there has been considerable interest in theories in which the SM fields are 
localized on a 3-brane in a higher dimensional spacetime. The idea that a multidimensional 
Universe may not be of KK type, but rather a low dimensional (3 + 1 in the case of our 
Universe) brane dates back to early eighties when independently Akama and Rubakov 
and Shapochnikov J7[ suggested models of our Universe as topological defects and later 
Visser || and Squires [f| described how particles can be gravitationally trapped on the 
brane. The possibility that gravity be trapped on a brane was suggested by Gogberashvili 


In papers of Antoniadis, Arkani-Hamed, Dimopoulos and Dvali |TT] , FJ, [13|] it was shown 
that the size of extra dimensions can be large, up to a mm scale. Then Randall and 
Sundrum proposed [[14] a scenario in the geometry is non-factorizable and even infinitely 
large fifth dimension leads to the Newton law at macroscopic distances along the brane. In 
this picture the concept of a massless graviton localized on a positive tension brane was of 
utmost importance. Another important element was a negative tension brane. 

Then in was suggested in [Q] and later and independently in that in some modification 
of RS scenario one can have even more strange behaviour - namely gravity is modified at 
both large and small scales. Even if explicit realizations of this idea suggested by Oxford 
group (KMPRS) |] and Gregory, Rubakov and Sibiryakov (GRS) p[ are different one can 
see that they actually related to each other and share several key elements - for example the 
presence of free negative tension branes. There is another common feature in both models 
- exponential suppression and tunneling effect (in case of GRS model it was discussed by 
Csaki, Erlich and Hollowood in 0). In this letter we demonstrate that these two models 

are actually two limiting cases of a more general model which we shall call H h or 

"multigravity" model and which interpolates between "bi-gravity" [p]] and "quasi-localized" 
gravity 0, 0. Besides this we are going to discuss the mechanism which will prevent 
negative brane destabilization - which is crucial for all these models. At the end we are 
going to discuss an important observation of Dvali, Gabadadze and Porrati Q that in this 
models four- dimensional graviton is massive and this leads to discrepancy ]I5| between 
Newton law and light bending (and other experimentally observed post-Newtonian effects, 
for example Mercury perihelion precession). The remarkable resolution of this potentially 
disastrous discrepancy is that in a full theory there is actually a smooth limit for graviton 
mass m — > if one take into account constraint / dy V—G^ (T£ — 2T|) = which is a 
necessary condition for stable compactification and correct cosmology on a brane || 

Let us first discuss the KMPRS model [p]] which consists of two positive branes located 
at the fixed points of a Si/Z 2 orbifold with one negative brane which can move freely in 
between (see Fig. [[]). It is easy to see that the two-brane RS model is nothing but the 
limiting case of our three-brane model, when the negative brane hits one of the positive 
branes. The model has two parameters, (at a fixed bulk curvature) - the warp factor and 
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Figure 1: H h model with two "+" branes at the fixed points and moving "— " branes. 

In the limiting case when i^Owe have a RS configuration. 



the x factor - which effectively measures the distance between one of the positive branes 
and the negative brane. The RS model corresponds to the limiting case x — 0. 

The GRS model can be obtained from this model by cutting negative brane in half, i.e. 
instead of one "—" brane with tension —A one can take two branes with negative tension 
—A/2 each and then move them apart, as on Fig. 2 If one then move the second "+" brane 
to infinity together with one of the " — 1/2" branes we shall get precisely the GRS picture 
- and the space between " — 1/2" brane is, of course, flat. Thus we have a new model - the 

H h model which depends on 3 parameters - II and Ir which are the differences between 

coordinates of left "+" and " — 1/2" branes and right "+" and " — 1/2" branes, and / 

which is a difference between coordinates of the two " — 1/2" branes. This model consists 
of four parallel 3-branes in an AdS$ space with cosmological constant A < 0. The 5-th 
dimension has the geometry of an orbifold and the branes are located at y = 0, y\ = II 

and y% = li + I and 2/3 = 2/2 + Ir, where y and y 2 are the orbifold fixed points^ (see 

Fig.2). Firstly we consider the branes having no matter on them in order to find a suitable 
vacuum solution. The action of this setup is: 

S= ( d 4 x f m dyV^G{-A + 2M 3 R} - £ f d A xV^ -G^ (1) 

where G^l is the induced metric on the branes and Vi their tensions. Here we included 
negative y and we are looking for solutions invariant with respect to Z 2 symmetry y — > —y. 



The notation is the same as in Ref. 14 



At this point we demand that our metric respects 4D Poincaree invariance. The metric 



3 The requirement that we have orbifold fixed points is not really necessary for our analysis, which is 
much more general 
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ansatz with this property is the following: 

ds 2 = e~ 2aiy) V ^dx^dx u + dy 2 (2) 

Here the "warp" function a(y) is essentially a conformal factor that rescales the 4D com- 
ponent of the metric. It satisfies the following differential equations: 

K) 2 = k 2 , ot'^-JfL-Siy-Li) (3) 

where k = sj 2 ±m 3 * s a measure °f the curvature of the bulk. The brane tensions are tuned 
to V = -A/k > 0, V 1 = V 2 = A/2k < 0, V 2 = -A/k > 0. It is convenient to introduce 3 
dimensionless parameters 

x L = Ml, x r = kl R , x- = kl— (4) 
It is easy to see that when X- = we have H h model with the hierarchy fac- 
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Figure 2: GRS model can be obtained from the left part of H h configuration. The 

whole configuration can be considered as two GRS models connected via flat space. 



tor w = e( XR ~ XL > and parameter x on Fig.l is xr. In the limit of infinite x_ we have 

two infinitely separated GRS systems. Then it becomes quite clear that H h model 

smoothly interpolates between this two models, spectrum that follows from the dimensional 
reduction. Then we have to find the spectrum of (linear) fluctuations of the metric: 

2 



ds 2 



M 3/2 V 



dx^dx v + dy 2 



We expand the field h lJiV {x 1 y) in graviton and KK states plane waves: 

oo 

h, u (x,y) = Y,h$(x)¥ n \y) 



(5) 



(6) 



n=0 
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where (d K d K — m 2 n ) hffi = and fix the gauge as d a h^l = h!^ = 0. The function $?( n \y) 
will obey a second order differential equation which after a change of variables reduces to 
an ordinary Schrodinger equation: 

{-\d 2 y + V(y)} ¥ n \y) = ^¥ n \y), ¥ n \y) = ¥ n \y)e^ 2 (7) 

where potential V(y) is determined by cr(y). Qualitatively it is 5-function potentials (attrac- 
tive for and repulsive for "—" branes) of different weight depending on brane tension and 
an extra smoothing term (due to the AdS geometry) that gives the potential a "volcano" 
form. 

An interesting characteristic of this potential is that it always gives rise to a (massless) 
zero mode which reflects the fact that Lorentz invariance is preserved in 4D spacetime. This 
mode is normalizable for finite x_ and becomes non-normalizable in a GSR limit of infinite 
x_. The interaction of the linearized gravitons to matter localized on a brane located at 
some y is given by 

Ant = M E * {n) (y)h$(x)T^(x) (8) 

with Tju, the energy momentum tensor of the SM Lagrangian and f(y) is some universal 
function. From this expression the Newton potential on a brane is given by 

(\fr( n h 2 (v) e ~ mr r ^ 2 (v)e~ mr 



„>„ M 3 r J M 3 



where p(m) is the spectral density and it is discrete for x_ = and becomes continuous in 
a GRS limit of infinite X-. We shall present full solution of this model for all X- elsewhere 
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Let us remind that one of the striking predictions of H h model is the fact that the 

first KK mode can be very light and strongly coupled compared to the rest of the KK 
states. This light mode can be so light that the corresponding wavelength can be by order 
of 1% of the observable size of the Universe while the second KK mode is in submillimeter 
region. Surprisingly enough this situation is not excluded experimentally and we called it 
"Bi-Gravity" - in all experimentally analyzed regions gravitational attraction is due to an 
exchange of two particles - the massless graviton and ultralight first KK mode. Only at 
scales larger than 10 26 cm will the first KK mode decouple leading to a smaller gravitational 
coupling beyond this length scale. The reason the anomalously light KK mode exists is due 
to the fact that with more that one "+" brane there will be a bound state on each of them 
when they have infinite separation. At finite distances there is a mixing between the two 
localized states. One superposition is the true ground state while the other configuration 
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has non-zero mass, but the gap may be very small - it is given by a tunneling factor. In 
the case x >> 1 we obtained a reliable approximation to its mass 



mi = 2ke- XL ~ XR (10) 

For w = 1 we get mi = 2ke~ 2x and the wave functions ^(O) — ^i(O) with exponential 
accuracy. The masses of the other KK states in the above region are found to depend in 
a different way on the parameter x. The mass of the second state and the spacing Am 
between the subsequent states have the form: 

m 2 ~ ke~ x , Am « eke~ x (11) 



where e is a number between 1 and 2. Equations fll0|) and (11) show that, for large x, the 
lightest KK mode splits off from the remaining tower. This leads to an exotic possibility in 
which the lightest KK mode is the dominant source of Newtonian gravity! 

Cavendish experiments and astronomical observations studying the motions of distant 
galaxies have put Newtonian gravity to test from submillimeter distances up to distances 
that correspond to 1% of the size of observable Universe, searching for violations of the 
weak equivalence principle and inverse square law. In the context of the graviton KK 
modes discussed above this constrains m < 10 _31 eV or m > 10~ 4 eV. Our exotic scheme 
corresponds to the choice mi m l(T 31 eV and m 2 > l(T 4 eV. In this case, for length scales 
less than 10 26 cm gravity is generated by the exchange of both the massless graviton and the 
first KK mode (C is some constant) 

U(r) - C ( * (0))2(0) ( 1 + i^m^l ) + (e-™ n ~ 2C ^° m 1 (12) 
f/(r) - C M 3 [r + (*(0))2( ) r ) +U[e M 3 r {U) 

so gravitational constant is Gn = 2(7 ^^/^ . According to this picture deviations from 
Newton's law will appear in the submillimeter regime m,2r < 1 as the Yukawa corrections of 
the second and higher KK states become important. Also the presence of the ultralight first 
KK state will give deviations from Newton's law as we probe cosmo logical scales mir > 1 
(of the order of the observable universe) with instead of G^. The phenomenological 

signature of this scenario is that gravitational interactions will appear to become weaker 
(in general case by the factor w [0) for distances larger than l/m\ ~ 10 26 cm! 

The idea behind GRS construction was totally different - they did not have normalized 
modes, but rather continuous spectrum, however they have a "resonance" which effec- 
tively was due to the fact that that the negative brane created a tunneling factor (negative 
brane acts as a repulsive potential) which effectively leads to a resonance in a wave function 
^m(O) of gravitons 

^(Q)= 2 ' r2/A + °(^ 4 ) (13) 
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Figure 3: Behaviour ^^(0) in a H h model (discrete spectrum). For a general H h 

model the discrete spectrum density increases with the increase of x_ and spectrum becomes 
continuous in GRS limit. 



, 2c 1 r dx ( Tr , 

U(r) = T ~ F - / -^-exp --x (14) 



From ([K|) one can see 



M 3 TrJ x 2 + l 

and we see that for r « 1/T we recover the Newton 1/r potential with 

Gn = W^f (15) 

and for r > 1/T potential is 1/r 2 and gravity is modified at large distances. 

But this is just the same as we have on our discrete case (see Fig. 3) - the reason is of 
course that T ~ e~ x . Thus we see that in both limits it the same tunneling factor e~ x made 
modification of gravity possible - in discrete case it is exponentially small splitting between 
masses of massles graviton and first KK mode, in continuous case it is exponentially small 
virtual resonance width. 

Of course one get the same behaviour in other models when positive branes are separated 
by negative branes. For example one may have a crystal [[17|] and in this case we have bands. 
Potential will be proportional to 

„ e~ mr 

U(r) ~ / dmpim) (16) 

J r 

where p(m) is band density and for narrow band we shall have multigravity again. All we 
need is that the width of the first band Tx is much smaller than the separation between 



first and second bands (more details will be given in [16 
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Figure 4: A periodic array of positive and negative branes is a model of Crystal Brane 
Universe W7\ . 



Let's make some comments about the stability of the brane configuration discussed here. 
It involves a negative tension brane located between positive tension branes. Necessarily the 
negative tension (three) brane is free to move (in this configuration it cannot be located at 
a fixed point) and the question naturally arises about its stability. Being negative tension 
it will be advantageous to increase its volume and this may happen both by changing 
its location through the effects of the warp factor or through the brane bending. It is 
instructive to consider the first effect on its own by determining the effect of moving the 
position of the brane a distance 5 from an initial position midway between the positive 
branes. In this case Q the warp factor at the negative brane decreases by the factor e _ft<5 , 
where k is the five dimensional curvature. Thus the volume of the negative brane decreases, 
implying the midpoint is a position of equilibrium for the "flat" brane. Consider now the 
effect of deforming the brane over some transverse size I which may reasonable be taken 
to be the horizon size. The mean displacement of the brane about its equilibrium position 
is proportional to R — RcosO ~ R8 2 = 16, where 6 = l/R -C 1. Thus the warp factor 
decreases by the factor e~ 2cKie where c is some constant. The deformed brane stretches 
by an amount «(1 + 6 2 /2) in each transverse dimension so the net change in the brane 
volume is (e~ CKie (l + # 2 /2)) 3 . Since the horizon size, I, will always be much larger than the 
Planck length rs k _1 ), this naive argument suggests that the warp factor is the dominant 
effect and that the bending of the brane is energetically disfavoured. If this is true the 
brane configurations discussed in this paper may be stable structures and the modified 
gravitational interactions discussed above may indeed be possible . 

Finally let us discuss problems raised in J|]. Contrary to massless four- dimensional 
graviton which has 2 polarizations the massive one has 5 possible polarizations. The tensor 



structure for graviton propagator is given by 

1/2 (g» a g"P + g va g^\ - tg^g^ 



c 



0(pp/p 2 



p 



(17) 

where for all non-zero m parameter t = 1/3 [|l5j but for m = t = 1/2. Then the one- 
graviton exchange amplitude for any two four-dimensional sources T* v and is given 
by 

T^G^Tlp = C^- 2 (18) 

where T = T£. Massive bodies at rest have only non-zero Too and T = Too- Newton 
constant can be found from amplitude 

C(l-t) 



~~ 2 2 00 1 00 

But there is another type of matter - light and for light T 
bending does not depend on t and is given by 

C 



Til rp2 



(19) 

0. The amplitude for light 

(20) 

The ratio of two amplitudes is 1 — t. Because both amplitudes are known experimentally 
it is easy to find t -and of course it must be t = 1/2 as in Einstein General relativity. If 
we take t — 1/3 the ratio of two amplitudes is going to be 2/3 instead of 1/2 - and this is 
factor 4/3 which was discussed in ||. 

However this is not the whole story. We have a five-dimensional General relativity and 
the graviton propagator in five dimensions has tensor structure HI 



G 



MN,PQ 



1/2 (g MP g NQ + g MQ g NP 



tg MN g PQ 



0(pp/p 2 



(21) 



where t = 1/3 i.e. as massive four-dimensional graviton. The reason is very simple - 
they both have 5 physical components. It is a well known fact in KK theory that after 
compactification the massless sector contains massless graviton, graviphotn and KK scalar 
- in total five degrees of freedom. At the same time massive KK excitations contain only 
massive 4d graviton. 

Let's then consider the full five-dimensional amplitude 

riMN,PQrp 
Tvr If 



T 



Mm 



PQ 



T 



MN 



1/2 (g MP g^ + g^g NP 



- tg MN g PQ 



T 



PQ 



22 



where T MN is a full stress-energy tensor. We don't have T^ 5 but besides T^ v there is T 55 . 
Then the total amplitude can be written as 

T^T^ - tTT + (1 - t)T|T 5 5 



2tTT 5 5 



T UV T» V - -TT 

5t-i5 



+ 



(1/2 - t)TT + (1 - t)T 5 5 T 5 & - 2tTT 



(23) 
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The first line is the correct four-dimensional amplitude - all known experimental data gives 
this amplitude, so it will be nice to see that the second line is identically zero. Amusingly 
at t — 1/3 the second line is a full square 

1 (T _ 2T 5 )2 (24) 

Why it is natural to have the condition (T — 2T|) = ? First of all this is written 
in momentum representation. All amplitudes under consideration have zero momentum 
transfer in the fifth direction. This means that we have to integrate over the fifth coordinate 
/ dyy/—Gffi (T — 2T|). But this is precisely the integral condition which was discussed in 
|H. The (755 remains in equilibrium and the size of an extra dimension stays fixed if the 
stress energy of the 5D matter satisfies the following constraint^]: 



(fx dy \f-G( 5 ) (T£ - 2T 5 5 ) = 0. (25) 

To show what does this condition mean let us consider standard Kaluza-Klein decom- 
position 

G&=f G r I), (26) 



e 2 -' 



and then consider the linear term 



SG { S N tMN = SGi uu + 5(e 27 ) T 55 (27) 



in any system. Using Gffi = e 1 g ilin one can rewrite the above expression as 

bGf 1N T MN = e-^g^ - 5 1 (T - 2Tf) (28) 

Thus, we see that is a source for the graviton while the combination (T — 2T|) is the 
source for the scalar (radion). In the presence of a mechanism that stabilizes the extra 
dimension (see for example J20| and |]2T| and references therein) the source for the constant 



mode of this scalar is identically zero. If we do not have a stabilization mechanism but, 
nevertheless, we look for solutions with a Newtonian limit, we still have to require the 
absence of the above term. 

The conclusion is that there is a smooth limit m — > in multigravity. The situation 
here similar to the Higgs mechanism. There is no smooth limit when the mass of a vector 
particle is given by hand - but there is a smooth limit when it emerges dynamically due to 
spontaneous symmetry breaking. 
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Note added in proof: During preparation of this paper for publications several 
papers were published in which the 4/3 problem [|J was discussed p2| , but not the constraint 
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